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Abstract. We show analytically that the perturbative expansion for the free energy 

2 ' of the zero dimensional (quenched) disordered Ising model is Borel-summable in a 

P^ , certain range of parameters, provided that the summation is carried out in two steps: 

n-^ ' first, in the strength of the original coupling of the Ising model and subsequently in 

r^ . the variance of the quenched disorder. This result is illustrated by some high-precision 

O I calculations of the free energy obtained by a straightforward numerical implementation 

I i ' of our sequential summation method. 
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1. Introduction 

One of the simplest cases to study tlie effects of impurities in magnetic systems occurs in 
some doped uniaxial antiferromagnets such as Feo.46Zno.54F2 ||T|, where the Zn impurities 
do not induce competing interactions and Monte Carlo simulations with the three- 
dimensional random-site Ising model reproduce accurately the behavior close to the 
Curie temperature f^ . Nevertheless, the successful apphcation of renormalization group 
techniques |^ to the calculation of critical exponents for the pure Ising model makes 
perturbative expansions a very appealing alternative to the Monte Carlo calculations 
mentioned above. 

The appropriate field theory for this problem is the usual u(f)'^ theory in three 
dimensions modified with a random mass term, and the analytic difficulties raised 
by the quenched impurities are usually bypassed with the replica trick: the system is 
temporarily replaced by n non-interacting copies for which the disorder is now annealed, 
yielding an 0{n) theory with coupling constant proportional to minus the variance of the 
quenched disorder w, with an additional cubic anisotropy term with coupling constant u. 
This theory can be studied perturbatively and the results corresponding to the original 
system recovered (in principle) in the limit n —>■ 0. 

Indeed, there are rather accurate analytic calculations of the critical exponents 
of the random-site Ising model [Q, but the situation is less satisfactory than for pure 
systems . For instance, the asymptotic parameter in disordered systems is y/e instead 
of e (where e = A — d, d being the dimensionality of the space), and the /3- function 
computed at three loops shows no stable fixed point |^. In an attempt to understand 
disordered systems in a simpler setting. Bray et al [H and McKane studied the 
asymptotic expansion for the free energy in the zero- dimensional case, with discouraging 
results: the toy-model was found non Borel-summable, in sharp contrast with the 
(ordered) zero- dimensional ucf)^ theory, where Borel-summability has been rigorously 
proved 0. (Indeed, Borel summability in three dimensions has been also proved [§.) 

We consider in this paper the same zero- dimensional problem of references P, |^, 
and show that a slightly more elaborated procedure can recover the free energy from the 
perturbative expansion: the Borel-summation has to be carried out first in the strength 
of the original coupling of the Ising model u and subsequently in the variance of the 
quenched disorder w. 

The layout of the rest of the paper is as follows. We first review the main results 
of Bray et al and McKane in section ^ the derivation of the Borel-summable double 
asymptotic expansion for the free energy is carried out in the third section; in section ^ 
we give the details of our numerical procedure and some illustrative examples, and the 
paper ends with a brief summary. 



z{i,,u)= r ^^-ii^+^)^'-i-^' (2) 
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2. Previous Results 

The free energy of the zero-dimensional Ising model with quenched dilution can be 
studied with a u(f)^ theory with a random mass term 

f{u, w) = - r -^ e-^'/^- log Zii^, u) (1) 

where 

dcj) 1- . -.■> 1 

is the partition function. Bray et al fix 

A = w/u > (3) 

and write the asymptotic expansion for the free energy as 

oo 

f{u,w)r^-J2AK{X)u'' (4) 

K=l 

where the ^^^(A) are polynomials in A with rational coefficients and degree K. Then 
they use the replica trick an a saddle point argument to infer that the asymptotic 
behavior of v4i^(A) as -ft' ^ oo is given by two terms: 

a^^\x) + aP{\), 0<A<1 

^ aP{X), 1 < A 

The term A)^ (X) is dominant for < A < | and Borel-summable, while 

A-(A) = - -^j^ exp(-7V^ + a) cos{fiVK + 6) (6) 

with coefficients 7, a, n and 6 that they ultimately adjust by comparison between this 
asymptotic formula and the numerical values of the Ak{X) for K up to 200 and several 
values of A greater than i. The increasing-period cosine oscillations suggest them that 
there is an essential singularity of the Borel-transform in the positive real u axis, and 
therefore the series is not Borel-summable for any finite disorder (although the term 
^K (X) is subdominant in the region < A < |). Despite the non-summability of the 
series @, i.e. acknowledging that there would be some error even if they could compute 
with all the terms, Bray et al set m = 1 and try to extract an accurate answer from @ 
using a conformal transformation suitable for oscillatory series, but they find that as 
soon as A > I the method fails. 

McKane [0] avoids the replica trick in favor of a more direct analysis that clearly 
identifies the two sources of non-analytic behavior of the free energy leading to the two 
contributions is equation (^: the branch cut along the negative u axis and the zeros 
of the partition function (0). Moreover, he obtains the exact values of 7, a, fi and 6 in 
terms of the zeros of the partition function and guesses the form of the nonperturbative 
contributions to the free energy that yield the asymptotic behavior (^. 

The possibility of alternative summation methods, however, remained unaddressed 
in these works. 



(7) 
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3. Summable asymptotic expansion for the free energy 

For analyticity considerations, we remark tliat in terms of the new variable 

_ 1 + V^ 

^ ~ 2mV2 
the partition function (j^) can be written as 

^W'.«)- ^1751^171 /. *e-"-J' (8) 

which in turn can be expressed in terms of a parabohc cyhnder function or a modified 
Bessel function: 

Z{^, u) = u~'/\2ny'/hh' z'/^Ki/^ilz^) (9) 

= (2«)-i/^e2^'D_i/2(v^z). (10) 

There is, however, another expression particularly suitable for later summability 
considerations: via a Mellin-Barnes integral [^, the partition function can be written 
as a Laplace transform in z"^, 

1 r°° dt e~^ * 



z r°° dt e 



-t 



{Rez>0). (12) 



(4u)l/4r(i) 7o [t(t + z2)]3/4 

Equations (||) or (|10|) show that 

F{z)=u^/^Z{ij,u) (13) 

is an entire function of z, and from equation (^ and reference []TU[ it follows that the 
only zeros of F{z) appear in complex-conjugate pairs Zk,Zk in the Kez < half-plane, 
and are given asymptotically by 

zl r^ -l\og 2 - i^{4k + 3) {lmzk>0,k = 0,1,2,...). (14) 

To stay close to the notation of Bray et al and McKane and yet show that we are dealing 
with a two-variable problem, we set 

A = w/u (15) 

/i = l/{2u"') (16) 

^{li,\) = f{u,w). (17) 
In terms of these variables, the free energy can be written as 

- $(/i. A) = ilog(2/i) + ^^-L_|_^dze-(^-'^)'/Mog[F(z)] (18) 

were the integral in the right-hand side is taken along the real axis, and is well defined 
for any complex yU and Re A > 0. Note that as z varies along the real axis in the 
increasing sense, F{z) traces the positive axis in the decreasing sense, but as we shift 
(say) upwards the integration path and z picks up a constant (positive) imaginary part 



Summability of the perturbative expansion for a disordered spin model 



CD 





II 


1 1 1 1 1 1 1 1 1 


1 


-Ill 


1 1 1 1 1 1 1 1 




' ^ 




- 




- 


- 






- 


0.5 


- 


z„ + 0.02 i 

Zo - O.OS i 


- 


- 


Z, 

2l 

^1 


+ 0.02 
- 0.02 


_ 




- 


;-■-;,'''' 


- 


- 


/ / 




- 


0.0 










/ 






;■ / 


s 

( 


" 


- ■ ;' 


\ 




- 


0.5 


- 


A- 


- 


— 


i 




- 



0.0 0.2 0.4 0.6 0.0 0.2 0.4 0.6 

Im F{z) 



Figure 1. Paths in the complex plane followed by the function F{z) as z varies along 
straight lines parallel to the real axis slightly below, on and slightly above the first two 
zeros of F(z) in the upper half plane (z — Zk + t, z — Zk He + t, — oo < t < oo). The 
dip in the right plot is a trace of the first loop of the spiral (shown in the left plot) 
which unfolds via a cusp (when F'{z) = 0) as the integration path shifts upwards. 



[i.e. z = t + ie,t E (—00, 00)], equation ([10|) shows that F{z) traces a spiral ending at 
the origin for t = 00. In particular, ii z = t + z^ the spiral passes again through the 
origin (see figure 1). Therefore, if we shift vertically the integration path to z = t + fi 
(with < Im/i 7^ Imzfc), we pick up a finite number of contributions from the zeros 
of F{z) whose evaluation in terms of the complementary error function is given by the 
following self-explanatory chain of equations: 



"1 pOO--{- u, 

$(/., A) = I log(2/i) + j—^ / dz e-(^-'^)VA iog[F(^)] 

+ 1: i^/_ d.e-<--'VMog|F(.)l 

= i log(2/.) + — — j^ dt e-* /' log[F(t + /i)] 



(19) 



(A 
E 

0<Imzj.<Im/^ 



27ri 



(A7r)i/2 



^^g-(.,-t-M)VA 



(20) 



1 1'°° 9 /, 

i log(2^) + — — J^ dt e-* /^ log[F(t + /i)] 



+ Itt ^ erfc 

0<lm Zf;<lm fj, 



—00 
fi 



Zk 



AV2 



(21) 



There are three key points to note in these expressions. First, the determination 
of the logarithms: they are real for real argument (as t -^ cxo) and must be followed by 
continuity according to figure 1. Second, for | Im/i| < Im^o there are no contributions 
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Table 1. Comparison among the K ~ 3 partial sum of equation (^, the shifted 
integral (P2[), the shifted integral plus corrections (pi|), and the unshifted integral ( [l8|) 
for A = 1 and three values of /i with < Im/i < Imzo, Imzo < Im/i < Imzi and 
Im zi < Im /i < Im Z2 respectively. 



^^ 



/i = 1 + 2i 



12=1 + 5i/2 



Partial sum 

Equation 
Equation (^1 
Equation flT 



0.0155640 -0.00752925 - i0.0099890 

0.0155379 -0.00752689 - i0.0099516 

0.0155379 -0.00620929 - i0.0107989 

0.0155379 -0.00620929 - i0.0107989 



-0.00624927 - i0.00593754 
-0.00624739 - i0.00593804 
-0.00908002 - i0.00573683 
-0.00908002 - i0.00573683 



from the zeros and we do have 

- $(/i, A) = I log(2/i) 



J —OQ 



(22) 



(A7r)i/2 J_, 

We stress that the right-hand side of equation ( ^21) is well defined whenever Im/i ^ Ira Zk, 
but it is equal to the free energy only for | Im/i| < Im^o- As a numerical illustration of 
this point, in the last three rows of table 1 we show numerical calculations of the shifted 
integral (0), the shifted integral plus corrections (plf), and the unshifted integral (0) 
for A = 1 and three values of fi with < Im/i < Irazo, ImzQ < Im/i < Imzi and 
Im Zi <Im fi < Imz2 respectively. Third, equation (p2| ) can be rewritten as a Laplace 
transform in A^^ 

-$(/i,A) = ilog(2/i) + ^y^ dte-*/^5(/i,t) (23) 

where 

1 



B{fi,t) 



log F{fx + t'/^)F{fx-t' 



/2^ 



(24) 



(47rt)V2 

The singularities oi B{fi, t) as a function of t are readily calculated in terms of the Zk and 
(apart from the trivial singularity at the origin) stay at a finite distance of the positive 
real axis whenever Im /i ^ Im Zk ■ 

Equations ( |12|) and ( |13|) imply that we can obtain a Borel-summable asymptotic 
expansion for F[z) and its derivatives, and equations (p3D and (|2^) that, as far as we 
stay away from the zeros of F, we can also obtain a Borel-summable series in A. Our 
final result will be, therefore, a Borel-summable series (in A) with each of its coefficients 
given by a Borel-summable series (in /t~^ = 4m) . 

Now we proceed to the details of the derivation. We first expand \og[F{z + /i)] 
in equation (p^) in (convergent) Taylor series around z = — or equivalently B{fi,t) in 
equation (p3|) around t = — and integrate term by term to get 

oo \n 

- $(/., A) ~ ilog(2/.) + E(log°^)^'"H/^)7^- (25) 



n=0 



'n! 



Then, we use the Borel-summable asymptotic expansion for F{z) that follows directly 
from equation ([12D 



Fiz) ~ {2z)-'/\Foii h 



v-2\ 



-f < arg^; < fj 



(26) 
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to obtain asymptotic expansions for the derivatives (logoF)(^")(/i). Incidentally, we 
mention that the summability sector of F{z) in terms of u is just — tt < argu < n, i.e. 
we are rederiving the well-known summability of the (ordered) zero- dimensional u(f)^ 
theory as explained, for example, in reference [Q. 
Note that for the zero-th derivative we have 



log[F(/x)] ~ - i log(2/i) + log [2Fo(i, |; ; -fi~')\ (27) 

oo (_-\\ku 

= -llog(2;.) + E^-^ (28) 

fc=i /^ 

and the logarithm cancels the first term in the right-hand side of equation (|22|). 

(Equation ( P5| ) is in fact the definition of the coefficients bk)- Higher derivatives have 

pure asymptotic power series that can be computed easily in terms of the zero-th 

derivative series 



(logoF)(2")(^)~ — 



/^ 



2n 



T{2n) 



+ E 
fe=i 



T{2n + 2k){-lfbk 
T{2n) /i2fc 



The final result is the double series 

oo /' 1 ^fcJ, oo -| / \ 



$(/i,A) 



v- v~l) bk , v^ J_ 
Z^ ,,2k "•" Z^ ^1 



fc=l 



/i 



n=l 



n! 



4/i2 



T{2n) ^T{2n + 2k) 



-l)% 



k=l 



T{2n) 
0) 



/i 



2fc 



(29) 



(30) 



whose structure in terms of the original variables is (cqo 

oo / oo \ 

-f{u,w)r^Y.[Y.^n,ku']w'^. (31) 

71=0 \fc=0 / 

This double series is formally equivalent to equations (BIO), (Bll) and (B13) in 
reference (the only difference is that we have transposed the indices in Cn^k), but 
with a definite ordering that cannot be altered without further considerations: to take 
advantage of the summability we should first Borel-sum the series in u (see the precise 
definition in the next section) 



c„(m) = Borel ^ c„,fcM'' 

k=0 



and then in w 

- f{u,w] 

instead of performing first the finite sum 

K 



Borel ^ Cn{u)w^' 

n=0 



(32) 



(33) 



(34) 



n=0 



followed by some kind of summation in u. Although this last reordering preserves the 
asymptotic nature of the series to the function defined by equation ( p2D (as we also 
illustrate in the first row of table 1), the numerical and analytic evidence presented by 
Bray et al shows that it spoils the summability of equation 
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4. Numerical algorithm and results 

Of the two standard methods to achieve the analytic continuation implicit in the Borel- 
summation (conformal mapping or Pade approximants) we have chosen the second 
because it does not require a precise knowledge of the singularities of the function, 
although some Pade approximants may introduce spurious singularities close to the 



integration path. Equation (^91) suggests that an appropriate numerical Borel-sum of 
the u series is 



U 



1 r^^ 



where P^^''^\t) is the [p, g]-Pade approximant for 

P+q „ fk 

Equation (|3^) has the additional advantage over other equivalent forms that the Pade 
approximants and their zeros (which we use later) need to be calculated only once. 
Furthermore, to avoid problems with the numerical integration we expand the Pade 
approximant as a polynomial plus partial fractions, i.e. assuming that all the poles are 
simple 

and the integration in equation (^) can be carried out in terms of complete and 
incomplete gamma functions evaluated at the q poles t^ of p[^'^l(t): 

p-q 



k=0 



+ ^5^ t |||e-'./.(_t,)--/^r(-2„ + i, -t,/u). (38) 



We have performed our calculations with the same 200 terms of the expansion (|28|) used 
in reference 0, using the arbitrary precision numerical capabilities of Matheinatica. 
Figure 2 shows some typical results of this first stage of our double-summation method. 
For four values of u and the lowest eleven coefficients Co(m), . . . , Cio(m) we plot a measure 
A of the number of correct digits to the right of the decimal point, 

A = - logio Icniuf''^''^ - c„(M)S^"^"^ed| (39) 

as a function of the order p of the diagonal Pade approximants P^'P^{t) in equation (|35|). 
As a general trend, the precision of the summed coefficients increases with the order of 
the approximant, although with smaller slope for higher values of u (note that we have 
used the same scales in the four plots). The irregularities in figure 2 are often due to 
ill-conditioned linear systems of equations in the calculation of the Pade approximants, 
and are usually corrected using higher precision or higher order approximants. We have 
also checked that neighboring off-diagonal Pade approximants give similar results. 
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Figure 2. Number of correct digits to the right of the decimal point in the Borel- 
summed coefficients co(u), . . . , cio(w) (from top to bottom a.t p — 10) as a function of 
the order of the diagonal Pade approximant in equation (p5|). 



The second stage, the summation in w, is equivalent to the n = summation of 
the first stage. Now the Borel-transformed coefficients c„(u)/r(n+ |) are precisely the 
coefficients of the convergent Taylor expansion of t^^'^B{fj., t), and to illustrate how well 
the Pade approximants continue the Taylor series beyond the radius of convergence, we 
plot in figure 3 the exact integrand -B(|, t) corresponding to m = 1 in equation (p3D, and 
its approximations t~^^'^P^'P^{t) for p = 10, 15 and 20. The radius of convergence of the 
Taylor series is marked in the plot by the vertical line at t = |zo " |P ~ 6.54. 

In figure 4 (again with u = 1 and for several values of A) we plot the number 
of correct digits to the right of the decimal point in the Borel-summed free energy as 
a function of the order of the diagonal Pade approximant. We have performed the 
calculations with all the c„(m) of the same precision, so that the numerical error is due 
only to this second-stage summation. The first feature worth mentioning is the steady 
increase in the precision of the summed free energy (with a few exceptions that are 
related again to the Pade approximants), although the convergence is slower for larger 
values of A. Note that the plateau in the summation of Bray et al for A = 0.2 ends 
approximately at fifty terms of the series, while at p = 40 (i.e. with 81 terms of the 
series) we have steady numerical convergence; A = 0.5 corresponds to the switch between 
asymptotic behaviors in equation (1), while in our reordering it is not a distinguished 
point; finally, for A = 0.8 Bray et al are already unable to get a good result, while our 
summation still gives four correct digits for A = 5. 
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Figure 3. Exact integrand B{^,t) corresponding to u = 1 in equation (23) and its 
approximations t^^/'^P^P-'P^{t) for p = 10, 15 and 20. The vertical line marks the radius 
of convergence of the Taylor series for B{^,t). 
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Figure 4. Number of correct digits to the right of the decimal point in the Borel- 
summed free energy /(I, A) as a function of the order of the diagonal Padc approximant 
used in the second stage of the double summation. 



5. Summary 



We have analyzed the Borel-summabihty of the perturbative expansion for the free 
energy in one of the simplest disordered systems, the zero dimensional quenched diluted 
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Ising model. This expansion is a double asymptotic series in the original coupling u of 
the Ising model and in the variance w of the quenched disorder, and we have shown 
that the expansion is in fact Borel-summable by a sequential method: first we sum in 
the coupling u and subsequently in the variance w. This summability results are valid 
at least for — vr < argw < tt and Rew > 0, with the understanding that if we want to 
recover the free energy and not the integral (P^ ) , additional erfc terms are needed in the 
domains specified in section |^. 

A straightforward numerical implementation of our double summation method has 
permitted us to obtain (with the same number of terms of the original series) accurate 
free energies for values of the parameters well beyond the region where Bray et al were 
already unable to get an estimate. 

We expect that the insight provided by these zero dimensional results will be 
helpful in the understanding of the more complicated higher-dimensional case, despite 
the significant qualitative differences motivated by Griffiths singularities |jll| in the 



generalized sense of Dotsenko [|T^]. In fact we feel that a clear understanding of the 
Griffiths phase in the field-theoretical setting will be needed before the problem will be 
finally solved. 

Finally, we would like to mention as open problems the study of the negative A 
region, and a more thorough understanding of the solution of this problem using replicas. 
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